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ABSTRACT 

We use deep, five band (100-500 jim) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFRfir, of 
galaxies in the Bullet cluster {z - 0.296), and a smaller background system {z - 0.35) in the same field. Herschel detects 23 Bullet cluster members 
with a total SFRfir = 144 + 14 M© yr~^ On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher 
SFRpiR (21 galaxies; 207 + 9 M© yr~^). SFRs extrapolated from 24 yum flux via recent templates (SFR24^m) agree well with SFRfir for ~60% of 
the cluster galaxies. In the remaining ~40%, SFR24^m underestimates SFRfir due to a significant excess in observed S 100/^^24 (rest frame 5" 75/5 ig) 
compared to templates of the same FIR luminosity. 

Key words. Galaxies: clusters: individual: Bullet cluster - Galaxies: star formation - Infrared: galaxies - Submillimeter: galaxies 
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1. Introduction 

In the last decade many studies have attempted to quantify the 
star formation rate (SFR) within cluster galaxies. Ultraviolet and 
optical observations have successfully identified trends between 
unobscured star formation and local environment, suggesting 
that star formati on in cluster c ore gala xies is genera lly more 
quenc hed (e.g. iKodama et al.l 120041: iPorter & Raychaudhuryl 
l20Q7h . However, star formation can be obscured by dust, which 
re-emits stellar light in the far-infrared (FIR), p eaking at a rest 
frame /Iq -100 ii\x\. Mid-infrared surveys (e.g. iMetcalfe et al.1 
120051: iGeach et a l. 2006; Fadda et al. 2008) have explored ob- 
scured star formation by estimating total FIR luminosity from 
template spectra. These templates are often based on small num - 
bers of well constrained local galaxies, e.g. lRieke et al.l (l20Q9l) . 

The PACS (Pogli tschetaD 12010") and SPIRE ("Griflan et alj 

2OI OI) instrun ients, onboa r d the ESA Herschel Space 
Observatory (iPilbratt et al.l |2Q1Q|) . enable unprecedented 
multi-band coverage of the FIR. The Herschel Lensing Survey 
(HLS; PI: E Egami) consists of 5-band observations (100-500 
yum) of 40 nearby clusters (z ~ 0.2-0.4). Nominally devised 
to exploit the gravitational lensing eflTe ct of massive clus ters 
to observe high redshift galaxies (see lEgami et aP |2Q1Q| for 
details on survey design), a useful by-product is deep FIR 
observations of the clusters themselves. At these redshifts, 
Herschel photometry spans the peak of the dust component. 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with 
important participation from NASA. Data presented in this paper 
were analyzed using "The Herschel Interactive Processing Environment 
(HIPE)," a joint development by the Herschel Science Ground Segment 
Consortium, consisting of ESA, the NASA Herschel Science Center, 
and the HIFI, PACS and SPIRE consortia. 



allowing an accurate constraint of far infrared luminosity, Lfir, 
and hence obscured SFR. 

During the Herschel Science Demonstration Phase, HLS 
observed the Bullet cluster (1E0657-56; z = 0.296). The rea- 
son for this choice was two-fold. First, previous studies re- 
port brig ht submillimeter galaxies in the background (e.g. 
iRex et al.l 12009.) . with HLS analysis presented in Rex et al.l 
(|2Q1Q|) . Second, the Bullet clust er is a recent collision of two 
clusters (JMarkevitch et al.l l2QQ2h . oflTering a unique laboratory 
for the study of star formation within a dynamic environment. 
The sub-cluster has conveniently fallen through the main clus- 
ter perpendicular to the line of sight (< 8° from the sky plane; 
iMarkevitch et al.ll2004l) . Analysis of X-ray emission shows that 
a supersonic bow shock precedes the hot gas, while the weak 
lensing mass profile indicates that this X-ray bright compo- 
nent lags behind t he su b -cluster galaxies due to ram pressure 
(IMarkevitch et al.1 l2QQ2l: iBarrena et al.1 120021) . A recent mid- 
infrared study by IChung et"aD (l2QQ9l) concluded that ram pres- 
sure from the merger event had no significant impact on the star 
formation rates of nearby galaxies. We can re-evaluate these pre- 
vious studies by using Herschel data to constrain Lfir directly. In 
this letter, we present an exploration of obscured star formation 
in this cluster environment. 



2. Observations 

2.1. Photometric data 

Five band Herschel imaging was obtained using two instru- 
ments: PACS (100,160 yum) covering approximately 8'x8' and 
SPIRE (250,350,500 yum) with a wider ~17'xl7' field. We also 
use Magellan IMACS optical, Spitzer IRAC and MIPS 24 pm 
maps with similar coverage to SPIRE, and high resolution HST 
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Fig. 1. Distribution of spectroscopic redshifts (0.27 < z < 0.37) 
for galaxies within the Bullet cluster field (outline). Herschel 
detected galaxies are also shown (filled). In addition to the Bullet 
cluster (z = 0.296), there is a background system at z = 0.350. 
Dotted lines show our membership limits of 3000 km s"^ and 
2000 km s"^ respectively. 



ACS images of the central 4'x4'. lEgami et aP (l2010h provides 
details of all the data, and presents Herschel FIR maps. 

The deep SPIRE maps have detection limits well below the 
instrument confusion limits. To avoid compiling sourcelists from 
confused maps, Herschel fluxes are measured at all Spitzer MIPS 
24 yum source positions. For a typical galaxy SED at z ~ 0.3, the 
24 yum map is much deeper than SPIRE, so even with a relatively 
high S/N > 10 cut (flux limit -100 yuJy), we can assume the 
inclusion of all sources contributing significant FIR flux. The 
use of mid-infrared source positions has the added advantage of 
decreasing the significance of flux boosting, which has not been 
addressed in this study. 

Photometric analysis followed the same procedure in all 5 
Herschel bands. An average PSF, measured from the brightest 
unblended sources in the image, was simultaneously fit to all 
positions in the 24 jim catalogue (without re-centering) using 
Daophot Allstar. At the longer SPIRE wavelengths, there is a 
higher probability of more than one 24 yum source falling within 
the FIR beam. In these instances, the objects are grouped to- 
gether at the 24 jim S/N-weighted mean position, treated as a 
single source, and flagged (see following su b- section). Fo r more 
details on the photometry technique see .Rex et al.l (l2010h . 

2.2. Spectroscopy and sample selection 

The spectroscopic redshift catalogue combines observations 
from three campai gns: Magellan IMACS multi-slit (856 targets; 
IChungetal.1 120101 Chung et al. in prep), CTIO Hydra multi- 
fiber (202; Fadda et al. in prep) an d VLT FORS multi-slit (14; J 
Richard, private communication). lEgami et al] (2010) provides 
further details. The merged catalogue comprises 929 sources 
within the SPIRE field. 

Figure [T] presents the distribution of spectroscopic redshifts 
for the range 0.27<z<0.37. An important aspect of this study 
is confidence in the cluster membership of galaxies. The Bullet 
cluster distribution peaks at z = 0.296, and we limit membership 
to ± 3000 km s"^ (0.286 <z< 0.306). In addition, this study also 
analyzes galaxies from a system at z = 0.350 in the same field, 
limiting membership to ± 2000 km s"^ (0.343 <z< 0.357). The 
systems have 362 and 95 known members respectively. 

The sample for this analysis consists of MIPS 24 yum sources 
with spectroscopically confirmed cluster redshifts. These two 
catalogues were merged by identifying the closest 24 yum source, 
within the RMS pointing error of MIPS (1.4'0, to the spectro- 
scopic position. For sample members grouped during the FIR 
photometry (previous sub- section), we examined the optical and 
IRAC colours of each group member, identifying the likely 



source of the mid- and far-IR flux. In cases where the sample 
member was not considered to be the source, or when the situa- 
tion was unclear, the object was rejected from the sample. 

In the final sample, there are 47 confirmed Bullet cluster 
members, and an additional 28 sources in the z = 0.35 system. 
Of these, 23 and 21 galaxies respectively are detected in the 
Herschel bands, highlighted by the filled distribution in Fig. [T] 
The background system has a much higher fraction of Herschel 
detections than the Bullet cluster (75% of 24 yum sources, com- 
pared to 50%). 

3. Results and discussion 

3.1. Far-infrared (FIR) spectral energy distributions 

For each source, the FIR spectral energy distribution (SED) is 
fit to all available Herschel data points, taking into account the 
upper limits for non-detections. The dust component is modeled 
by a modified, single-temperature, blackbody 



Sy = Niylv^fBAT) 



(1) 



where 5'y is flux density, yS is dust emissivity index (fixed at 1.5; 
using yS = 2.0 would vary Lfir by < 15% on average) and By{T) 
is the Planck blackbody radiation function for a source at tem- 
perature T. The shape of this optically thin (rather than thick) 
blackbody imitates the inclusion of a secondary (warm) dust 
component. As we are concerned only with Lfir and SFR, the 
parameterization of the data is the most important aspect, and T 
is used purely as a fit parameter. Galaxies within the PACS field 
have well constrained fits, and T is allowed to float freely. For 
those without PACS data (-40%), To has been forced to a narrow 
range centered on the mean value from the constrained SEDs (30 
± 1 K). Forcing To to a similarly narrow range about values ~1 cr 
from the constrained mean, varies Lfir by < 25%. Bias in Lfir 
due to model priors is comparable in scale to systematics from 
instrument calibration. 

Lfir is integrated over (rest f rame) X^ = 8-10 00 yum from 
which SFRfir is derived using the lKennicuttI (11998 ) relation. As 
an illustration. Fig. [2l displays the FIR SED fits for five of the 
most luminous galaxies in the sample. Thes e simple fits may un- 
derestimate Lfir by up to a factor of 1.8 (iRex et al .1120101) . as 
they lack a mid-infrared component. Future analysis will fully 
account for additional components. For the purposes of this 
study, a blackbody fit is sufficient. The luminosities of galax- 
ies in the background system have been de-magnified using the 
Bullet cluster lensing model of Paraficz et al. (in prep). The re- 
maining figures in this paper present de-magnified values. 



3.2. Star formation rates within the systems 

The system at z = 0.35 contains IR galaxies brighter than those 
in the Bullet cluster, with three galaxies meeting the LIRG cri- 
terion [log(LFiR/Lo) ^ 11.0] and an additional two within 1 cr. 
A further 10 members have log(LFiR/LQ) > 10.5. In contrast, the 
Bullet cluster contains two LIRGs, and only six other galaxies 
brighter than log(LFiR/Lo) = 10.5. 

The total star formation rate of the 23 Bullet cluster galaxies 
is 144 ± 14 Mq yr"^ The 21 galaxies in the background system 
are, on average, -50% more active, with a total SFR = 207 ± 
9 Mq yr"^. Only five of these galaxies have been magnified by 
more than 20%, and the minimum detected SFRs are similar in 
each system. Therefore, it is unlikely that the higher total SFR 
is due to a decreased lower limit caused by magnification. The 
diff'erence is likely to reflect the mass of the systems, although 
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z = 0,297 
L^iP = 8.50e+10 L, 
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Fig. 2. Photometric data for five of the most FIR luminous galaxies in the sample. Blue = R\CS; red = SPIRE; grey = IRAC/MIPS. 
Flux densities are as observed (i.e. not de-magnified). Redshift and, for background system galaxies, magnification factor, yu, are 
displayed at the top of each panel. Lfir and SFRfir derived from the best fit blackbody (black line) are also shown and have been 
de-magnified where necessary. 
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Fig. 3. Smoothed density maps for IMACS 5-band flux (left pan- 
els), 24 yum flux (central) and SFR calculated from Herschel data 
(right). The sources are binned by confirmed system member- 
ship, with the number of contributing galaxies displayed in the 
upper-right of each panel: upper row for Bullet cluster (over-plot 
in orange by the weak lens mass map); lower row for z = 0.35 
system. All maps are Gaussian smoothed to the SPIRE 250 jim 
beam size (18'' FWHM). 



the lower SFR in the Bullet cluster may indicate that cluster- 
cluster mergers are not important for triggering FIR starbursts. 

Figure [3] displays the spatial distribution of the Herschel- 
derived SFR for the two systems. Flux densities in optical B- 
band and 24 yum are shown for comparison. An initial exam- 
ination suggests that the Bullet cluster exhibits a radial trend 
in SFRfir (lacking significant FIR detection towards the cen- 
tre), reminis cent of that found in other contemporary studies 
(iBrag fia et al. 2010; iPereiraet al.1 12010) . The gradient in the 
Bullet cluster SFR is examined in detail in Chung et al. (in prep). 

The IR and optical flux of the background system trace sim- 
ilar distributions, whereas in the Bullet cluster, the 5-band flux 
is more centrally concentrated, away from the IR sources. This 
may indicate a difl'erent trend in dust retention for the two sys- 
tems. While the 24 yum and Herschel SFR density maps generally 
trace the same distribution, there are significant outliers: bright 
24 yum sources with relatively lower SFRs, and vice versa. In the 
following se ction, we compar e the SFR estimated from 24 yum 
(through the lRieke et"aDl2009l templates) to the SFRfir. 



3.3. 24/dm as a Lfir predictor in nearby clusters 

The mid-infrared bands, e.g. MIPS 24 yum, are often used to esti- 
mate far infrared luminosity, Lfir, and hence obscured SFR, via 
template FIR SEDs such as Rieke et al. (2009). Those authors 
provide a simple formula (their equation 14) to convert 24 jjm 
flux directly to SFR. The templates are based on local (U)LIRGs 
(z ^ 0.1), and at high redshift, may not be valid. Here, we test 
the template accuracy for cluster galaxies at z = 0.3. 

In Sect. 3.2 (Fig. [3]), we suggested that while 24 yum flux 
and SFRfir follow the same general distribution, they are not 
perfectly correlated. A direct comparison of SFRfir to SFR24^m 
(Fig. S plotted against the dust-peak-mid-IR flux ratio) leads 
to the same conclusion. For -60% of galaxies, the two SFRs 
agree well. However, there are several galaxies (-30%) that have 
severely underestimated SFR24^m, and these also display system- 
atically redder S 100/^24. If SFRfir is underestimated by the sim- 
ple blackbody fits (Sect. 3.1), the SFR24^m predictions are corre- 
spondingly worse. 

Are the under-predicted SFR24^m caus ed by t he redder 
^100/^24 colours? Figure [5] examines the Rie ke et all templates 
more closely, comparing them to the Herschel fluxes. For tem- 
plates spanning the Lfir range of the observations, the agree- 
ment is good for Aq ^ 200 yum. However, at 100 yum there are 8 
significant outliers; we define 100 yum excess galaxies as those 
with 5" 100/5' 24 (rest frame Sjs/Sis) > 30 as the templates pre- 
dict S 100/^24 ~ 20. Only templates with very high luminosities, 
i.e. log(LFiR/LQ) > 12, match the observed S100/S24, but even 
the brightest sample galaxy has only log(LFiR/LQ) ~ 1 1.5, while 
most are log(LFiR/Lo) < 11. Although high Lfir templates have 
5 100/^24 ~ 30, their lower peak wavelength leads to an under- 
prediction at A ^ 200 yum in at least three observed SEDs. We 
also compare to the least active iDale & Heloul (l2002h FIR tem- 
plates (a = 1.8-2.5). Th e locus of these are substantially similar 
to the low luminosity Rieke et al. templates and thus also only 
under-predict5' 100/^24. We stress that, unlike Lfir and SFRfir, 
the presence of a 100 yum excess is independent of the blackbody 
fits and the systematic uncertainties therein. 

Galaxies with a 100 yum excess account for -40% of clus- 
ter members detected with PACS, and cover the entire range of 
Lfir sampled. Above a nominal luminosity limit of lO^^L©, 55% 
of Bullet cluster galaxies have the 100 yum excess. The fraction 
in the background system is lower at 36%. This may indicate a 
trend with environment, or could be due to the off'-centre view 
of the latter system (i.e a potential radial trend). High resolution 
HST imaging covers five of the eight 100 yum excess galaxies 
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Fig. 4. Ratio of S FRftr (from blackbody fit) to SFR24/im (via 
iRieke etaP 12009) versus the flux ratio 5'ioo/5'24. For galaxies 
outside the R\CS field, 100 yum is predicted from the blackbody 
fit. Dashed line is equality and shaded region indicates 50% dif- 
ference in the SFRs. All galaxies with under-predicted SFR24;^m 
have redder S ioo/5'24. 
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Fig. 5. Observed Herschel fluxes normalized by 24 yum for galax- 
ies in b oth systems. Symbols as described at top-left. Rieke et al] 
(l2009h average templates for 9.75 < log(LFiR) < 11.0 (red- 
pink) plus one example high-LpTR templ ate (cyan). Locus of low- 
activity templates from iDale & Heloul (2002) {a = 1.8-2.5) is 
shaded green. All templates are normalized at /lobs = 24 yum (Aq 
= 18 yum). Templates in Lfir range of the observations under- 
predict S 100/^24 for 40% of sources. High Lfir templates do not 
match the shape at /I > 200 yum. 






Fig. 6. HST ACS thumbnails of five 100 yum excess galaxies 
(5" 100/^ 24 > 30; increasing from left) do not suggest a single 
source population. 



clei (AGN) which are not considered in the templates. Using a 
simple power law to parameterize flux in the range 24-100 yum, 
we estimate the AGN contribution to to tal bolometric lumino sity 
via the ^60/^25 indicator for ULIRGs (IVeilleux et al.ll2009L fig. 
36). None of the 100 yum excess galaxies have predicted AGN 
fractions > 30%. However, we may be under-predicting the con- 
tribution if the mid-IR SED steepens beyond 60 yum, or if the 
indicator breaks down for galaxies in this luminosity range. 

Herschel R\CS observations of z ~ 0.2 LoCuSS clusters 
(without the advantage of complementary SP IRE data), display 
a similar fraction o f 100 yum excess galaxies (ISmith et al.ll2010l: 
iPereira et al.l2010h . However, the high redshift fie ld sample from 
the HLS Bullet cluster observations (JRex et al.ll2010 ) lacks a 
comparable excess at Aq ^ 75 yum. These results suggest that the 
eff'ect could be either redshift dependent or cluster- specific. HLS 
is well placed for further analysis of the S 100/^24 phenomenon, 
as the combined PACS -h SPIRE data ensures that both the excess 
and entire FIR component can be constrained simultaneously. 



4. Conclusions 

Using deep Herschel observations (100-500 yum) to fully con- 
strain the FIR component, we derive obscured SFRs for galax- 
ies in the Bullet cluster (z = 0.296), and a background system 
(z = 0.35) in the same field. Herschel detects 23 Bullet clus- 
ter members, with a total SFRfir = 144 ± 14 Mq yr"\ while 
the background system contains 21 detections but -50% higher 
SFR (207 ± 9 Mo yr"^). The relative distributions of SFRfir 
and optical flux suggest a difl'erence in dust retention between 
the two systems. For -60% of galaxies, SFRfir agrees well with 
estimated SFRs from 24 yum flux via recent templates. However, 
the remaining galaxies display a significant excess at 100 yum 
(Aq ^ 75 yum) compared to templates, which causes an under- 
prediction in SFR24;^m. We note t hat such an exces s is not found 
in the high redshift, field sample (iRex et al .1120101) . Future stud- 
ies will exploit the full range of 5 -band Herschel cluster observa- 
tions available in HLS, to form a more complete understanding 
of the environmental eff'ect on obscured star formation rates, and 
explore the origin and dependencies of the 100 yum excess. 
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ticipation by NASA. Support for this work was provided by NASA through an 
award issued by JPL/Caltech. 



(Fig. O. Despite the small number, the galaxies span a broad 
range of types and morphologies. Further examples are required 
for a firm conclusion, but these suggest that the 100 yum excess 
is not due to a single population of galaxies. 

The 5" 100/5' 24 colours alone may have led to the conclusion 
that the 100 yum excess was due to galaxies with generally colder 
dust. However, fits to the combined HLS R\CS-hSPIRE photom- 
etry suggest that this is not the case. Rather, the excess may be 
due to an additional warm dust component or active galactic nu- 
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